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Abstract 

We present a new analysis of the E772 and E866 experiments on the nuclear de- 
pendence of Drell-Yan (DY) lepton pair production resulting from the bombardment 
of 2 H, Be, C, Ca, Fe, and W targets by 800 GeV/c protons at Fermilab. We em- 
ploy a light-cone formulation of the DY reaction in the rest frame of the nucleus, 
where the dimuons detected at small values of Bjorken 2:2 < 1 may be considered 
to originate from the decay of a heavy photon radiated from an incident quark in 
a bremsstrahlung process. We infer the energy loss of the quark by examining the 
suppression of the nuclear-dependent DY ratios seen as a function of projectile mo- 
mentum fraction x\ and dimuon mass M. Shadowing, which also leads to nuclear 
suppression of dimuons, is calculated within the same approach employing the results 
of phenomenological fits to deep inelastic scattering data from HERA. The analysis 
yields —dE/dz = 2.73 ± 0.37 ± 0.5 GeV/fm for the rate of quark energy loss per unit 
path length, a value consistent with theoretical expectations including the effects of 
the inelastic interaction of the incident proton at the surface of the nucleus. This is 
the first observation of a nonzero energy loss effect in such experiments. 



1 Introduction 



Nuclei can serve as a unique tool to study the space-time development of the strong inter- 
action during its early stages, which is inaccessible in collisions between individual hadrons. 
The Drell-Yan (DY) reaction |J on nuclear targets provides, in particular, the possibility 
of probing the propagation of partons through nuclear matter in its ground state, with the 
produced lepton pair carrying away the desired information on the energy and transverse 
momentum of the parent projectile quark after it has traveled in the nucleus. In this pa- 
per, we are specifically interested in the determination of the rate of energy loss per unit 
length k = —dE/dz of a fast quark propagating through a nucleus. We believe that the DY 
reaction data provides the cleanest way to single out such energy-loss effects. 

One might also consider using other reactions for identifying energy-loss effects, with 
deep-inelastic scattering (DIS) coming to mind as one possibility. However, experimentally 
|fj, [3], |, H , it has proven difficult to identify the partonic energy loss, not to mention to specify 
how to best measure this potentially important phenomenon. The origin of the difficulty can 
be appreciated by considering the following experiment. Suppose one were to measure the 
nuclear (A) dependence of semi-inclusive DIS, with known 4- momentum transfer. In this 
experiment, the A dependence of the energy (momentum) carried by a hadron created by the 
struck quark is supposed to convey the information regarding parton energy loss. However, 
besides the vacuum (initiated by the DIS) and induced (due to multiple interactions of the 
struck parton) energy loss, one must also account for the hadronization parton — > hadron. 
As a rule, this hadron is subject to absorption when it is produced inside the nucleus. The 
complicated dynamics of this (see || and the comparisons of the predictions to data f5|) 
makes it extremely difficult to single out the net effect of energy loss. 

High-pr production of hadrons off nuclei is even more complicated, since it includes a 
convolution of the parton distribution in the incident hadron with the high-p^ parton-parton 
scattering cross section. Multiple parton interactions inside the nuclear medium (Cronin 
effect) also makes the interpretation extremely complicated. Yet another reaction sensitive 
to energy loss is charmonium production off nuclei. The violation of x<i scaling observed in 
the E772 experiment [[jj has already suggested the presence of energy-loss effects. However, 
final state absorption and coherence effects (see |J) do not allow a clear identification of 
energy loss in the data. 

In the DY reaction, the ratios Ra/d of the cross section for a heavy nucleus A compared 
to a light one (the deuteron D, say), are particularly relevant for the study of parton energy 
loss in nuclei. This is especially true in the region of large values of the the longitudinal 
momentum fraction of the projectile hadron carried by the DY pair, where significant nuclear 
suppression is clearly seen in this ratio. In an early study , it was shown that a ratio such 
as Ra/d would be sensitive to k because the fraction x q of the light-cone momentum of the 
incoming hadron h carried by projectile quark q shifts to smaller values in hA collisions as 
a result of initial-state energy loss. In the case of the DY reaction, this shift Ax q suppresses 
the ratio Ra/d because it results in the sampling of the parton distribution at larger x q 
(for the same momentum of the dilepton), where the projectile parton densities are smaller. 
Thus, the ratio Ra/d must decrease towards large Feynman xf — > 1 where energy loss leads 
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to a strongest suppression. 

At the same time, nuclear shadowing was observed in Ref. |§ to be a competing source 
of suppression at large Xp (more specifically, small values of the momentum fraction x<i 
of the target). As it is easy to mix up energy loss and shadowing, one must take special 
precautions to disentangle these two effects when analyzing experimental data. 

The first Drell-Yan data suitable for such an analysis was obtained in Fermilab E772/E866 



experiments. An analysis of the E772 data |T0| was made in [|IT| ignoring shadowing and 
assuming that k rises linearly with energy. The latter assumption was criticized [12] for 
violating the Landau-Pomeranchuk principle. 

A better analysis was performed recently by the the E866 collaboration using the E866 
data |R|. The analysis of Vasiliev, et ai, which attempted to improve on that of |lTJ by 
including shadowing, considered three scenarios for energy loss, one of which was the same 
All three of these gave k consistent with essentially zero energy loss, in contrast 

Thus, shadowing was found in 



as m 
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to the value k m l.hGev / fm found in JTT 



13 to be the 



main source of nuclear suppression of the DY cross section at large xp. 

Our present work differs from the previous analyses in that we attempt to unambiguously 
separate shadowing and energy loss by calculating the shadowing correction to the DY data 
using theory. The concept of coherence length, which plays a key role in respecting the 
Landau-Pomeranchuk principle, is essential to this formulation. It is our belief that the 
procedure employed by Vasiliev, et al. overestimated the shadowing contribution and hence 
substantially underestimated k. The reason is that the shadowing correction was taken 
from a phenomenological analysis |L4| that had already attributed the suppression of Ra/d 
observed in E772 data at large xp entirely to shadowing. Preliminary results of our analysis 
are reported in Ref. |15| . 

We describe the interplay between quark energy loss and shadowing by working in the 
target rest frame using the light-cone approach of Refs. fL6 , |17| , |I8| |P9| , where these are 
given as separate contributions to the DY cross section. The target rest frame formulation 
is discussed in Sect. 2.1, where we relate the DY process to projectile fluctuations containing 
the DY pair. In terms of such fluctuations, the DY reaction may be viewed as occurring 
when interactions with the target break the coherence of the fluctuation and free a DY 
pair, i.e. bringing it on its mass shell. Such an interpretation is quite different from the 
conventional partonic treatment of DY process as annihilation qq — > 11. However, the 
partonic interpretation is known to be Lorentz noninvariant and vary from frame to frame. 

Such a rest frame formulation bears a close analogy to the more familiar light-cone 
description of deep-inelastic scattering (DIS) at small x pOp . In that case, the incident 
virtual photon develops a quark- ant iquark fluctuation which is brought to its mass shell by 
the interaction with the target. Again, this picture looks very different from the conventional 
parton model interpretation in the Breit frame, where it is described as absorption of the 
incident virtual photon by a quark (antiquark), which in this case belongs to the target. 

The two processes, DY and DIS are known to be closely related in the parton model via 
QCD factorization. It is not surprising that the rest frame descriptions of these reactions 
also look similar. One of the advantages of this approach is the clear and simple treatment 
of nuclear shadowing, which is described in terms of the usual Glauber-Gribov theory. 
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Shadowing in both processes, DY and DIS, is controlled by the same universal color-dipole 
cross section f2l| , which is the cross section for freeing the fluctuations. The strategy used 
in present paper is to treat the dipole cross section phenomenologically by fitting it to DIS 
data and then verifying the theory by comparing it to DIS scattering data on nuclei in 
the shadowing region (where there is no danger to mix up shadowing with energy loss). 
Then one can safely predict shadowing for the DY process. Such a strategy is based on the 
universality of the dipole cross section, which is a manifestation of QCD factorization. 

The lifetime of the fluctuations, or coherence length, is another crucial quantity for 
understanding how energy loss and shadowing occur in the target rest frame. The coherence 
length is discussed and calculated in Sect. [2.2| . This turns out to be mainly a function of 
Bjorken x 2 , with some corrections that dramatically deviate from QCD factorization towards 
the smallest x 2 . 

We discuss the use of the DY reaction as a probe for quark energy loss in Sect. |[ Section 
|3.1| gives an overview of model expectations for energy loss. The string model (Sect. |3.1.1| ) 
predicts a constant rate of energy loss following the first inelastic collision of the incident 
hadron on the nuclear surface. The magnitude of this is about the same as that of the string 
tension, k s « 1 GeV / fm. Multiple interactions of the quark in the nuclear medium lead to 
an additional induced energy loss whose rate rises linearly with the path length of the quark. 
Numerically this is a small correction to the dominant constant term. Similar effects follow 
from perturbative QCD (Sect. p. 1.2 ). The first inelastic interaction of the incident proton 
in the nuclear surface initiates a long-lasting gluon bremsstrahlung, providing a constant 
rate of energy loss of about the same value as that given by the string model. Energy loss 
induced by quark rescattering also rises quadratically with the length of the path and is of 
similar value as in the string model. These two sources of energy loss are complementary 
and must be added, giving an expected rate of energy loss of about 2 GeV j fm. 

In the target rest frame, the suppression of Ra/d by K is associated with short-lived fluc- 
tuations in which the DY pair is released immediately after the interaction of the projectile 
quark with a target nucleon. The produced dilepton pair carries undisturbed information 
about the energy of the quark traversing the nucleus only for such fluctuations. An impor- 
tant, although simple, consideration that arises here is the determination of the length of the 
path in nuclear matter over which the parton in the initial state loses energy. It is demon- 
strated in Sect. O that the mean length of this path does not follow an A 1 / 3 dependence 
but is, in fact, quite short compared to the nuclear radius. The distribution over path length 
is calculated for different nuclei. Nuclear suppression caused by energy loss is then given as 
a convolution of the quark distribution function of the incident hadron, modified by energy 
loss, with the DY cross section for the quark-nucleon interaction (Sect. |3.3| ). The former 
is borrowed from phenomenological models successfully describing soft hadronic collisions 
(Sect. |3.3.1|) , while the latter is fitted to data from the E772 experiment on a deuterium 
target (Sect. [3T3T2D . 

The other consideration in understanding the suppression of Ra/d is, as we have re- 
marked, the shadowing process. The critical quantities for describing shadowing in this 
approach are the coherence length and the effective fluctuation-freeing cross section (aver- 
aged over different fluctuations). The effective cross section is calculated in Sect. f| in terms 
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of a color-dipole cross section that describes data from HERA for the proton structure func- 
tion at high Q 2 . The shadowing correction is then calculated in an eikonal description of 
the multiple scattering of the incident quark in terms of the color-dipole cross section. In 
Sect. f| we also confirm explicitly our description of shadowing by comparing the theory to 
deep inelastic scattering data, where energy loss is not an issue. We make our determination 
of k from R*ju D (x\, M 2 ) based on the theoretical description of the effects of shadowing in 
Sect. |5]. We fit the parameter k to the entire set experimental data without preselection, 
and the result is in agreement with the expected value of k. 

We also performed a variety of tests to check the stability of our results. These tests 
are described in Sect. 5.2. In one, we examined the sensitivity of the rate of energy loss 
to the relative contribution of shadowed and nonshadowed events. We determined that 
selecting only those events with small x%, where shadowing is the dominant effect, does not 
change the results of the fit. We also artificially enlarged or eliminated shadowing, corrected 
shadowing calculations for multiple interactions, etc., and found stable results. These tests 
also provide a scale for the systematic uncertainties in our analysis. 

Section 5.3 is devoted to the important issue of the disagreement between our determi- 
nation of the rate of energy loss and that of a previous analysis of the E866 data for the 
DY reaction flT3"| , which detected no energy-loss effect. In Sect. |5.3.1| we discuss why the 



choice of the target rest frame is the natural one for formulating the theory of energy loss. 
In Sect. 5.3b we find that the different conclusions reached in Ref. [l3| can be explained 
by the substantially different space-time variation in the DY reaction at small 22 and the 
larger density of antiquarks at large x<i employed in this work. We justify our own results 
based on detailed theoretical arguments. 

In Sect. 6 we summarize the results of the present analysis. 



2 Drell-Yan reaction in the nuclear rest frame 

2.1 qq — > II annihilation or q — > qll bremsstrahlung? 

In the target rest frame, the Drell-Yan reaction at small X2 corresponds to the electromag- 
netic radiation of a lepton pair by a projectile quark (valence or sea, depending on the value 
of £1), rather than to qq — > 11 annihilation |16| |. In this case, the // pair is imagined to be 



liberated from a virtual fluctuation of the projectile when it interacts with a nucleon of the 
nucleus. Two examples of Feynman diagrams contributing to this are shown in Fig. [1|. The 
cross section corresponding to these and other diagrams not shown has a factorized form in 
the impact parameter representation [JH], [T7|, [T^| . 

t^ Y d = J d2rT \^gii( r T, a)\ 2 a qq (ar T , s qN ) , (1) 

where s q N is the square of the center-of-mass energy of the quark-nucleon collision, M is 
the dilepton effective mass, and ^ q u(r T , a) is the light-cone wave function of the \qll) Fock 
component of the projectile quark. This wave function depends on both the transverse 
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Figure 1: Feynman diagrams for the cross section of lepton pair produc- 
tion in a quark-nucleon collision with the 7* radiated before (left) and 
after (right) the interaction with the target. 



separation tt between the 11 and the recoil quark and the fraction a of the initial light-cone 
momentum carried by the // pair. Here, <Jq q (p, s) is the universal dipole cross section for the 
interaction of a colorless dipole qq with transverse separation p = arx from the nucleon, 
introduced in | 2l| . The appearance of the color dipole cross section in Eq. ([!]) arises from 



the difference between the interaction for \q) and \qll) fluctuations. This means, among 
other things, that the DY cross section receives no contribution if the transverse positions 
of the initial and final quarks coincide, and at small p this cross section vanishes oc p 2 . 
Said otherwise, when — > 0, the strong interaction cannot discriminate between the Fock 
components \q) and \qll) and therefore radiation is not possible. 

As usual for bremsstrahlung, one cannot say whether the (virtual) photon is radiated 
before or after the interaction with the target (via gluon exchange); both are represented in 
Fig. |]. This quantum-mechanical uncertainty in the time of radiation is called the coherence 
time, which is the same as the coherence length since we are near the light cone. In terms 
of the light-cone approach, the coherence time can also be interpreted as the lifetime of the 
\qll) fluctuation of the incoming quark. 

As already noted, one may distinguish this approach from the usual partonic interpre- 
tation of the DY reaction noting that the space-time development of high-energy reactions 
treated in terms of the parton model is not Lorentz invariant and depends on the reference 
frame. Indeed, since even the time ordering varies from frame to frame, the lepton pair 
produced at small values of Bjorken 22 "C 1 in the DY reaction may be equivalently viewed 
as annihilation qq — > 11 in the rest frame of the photon |IJ or as the decay of a heavy photon 
radiated in a bremsstrahlung process. 
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2.2 Coherence length 

The coherence length is an important quantity controlling nuclear effects in the DY reaction. 
If the coherence length is longer than the mean internucleon separation in the nucleus, 
different bound nucleons compete in freeing the 11 fluctuation. This phenomenon is known as 
nuclear shadowing. On the other hand, if the coherence length is very short, the fluctuation 
has time to interact with only one bound nucleon. In this case, all nucleons contribute 
equally to the DY cross section, i.e. there is no shadowing. 

The coherence length for a fluctuation of a projectile quark q — > 11 q is given by the 
energy denominator, 

c " Ml -ml' 1 } 

llq q 



where 

Ml = 

1 — a a ail — a 



r2 _ M 2 m\ k\ 



M i q = + — + 7777^3 ( 3 ) 



is the effective mass squared of the llq fluctuation; p+ = E q +p^ is the light-cone momentum 

of the incident quark; M is the dilepton effective mass; and, kr and a are the transverse 
momentum and fraction of the light-cone momentum of the parent quark carried by the 
pair, respectively. 

If energy is conserved, the longitudinal momentum transfer between the initial state \q) 
and the fluctuation \qll) is q c = l/l c . Thus, one can say that the coherence length is the 
maximal longitudinal distance between fluctuations that are in phase. 

The DY variables Xi and x 2 satisfy the equations, 

x 1 x 2 = — 1 (4) 

X\ — X2 = Xf , (5) 

where xf is the Feynman variable, and they are interpreted in the parton model as the 
Bjorken variables of the annihilating quark and antiquark. They can equivalently be defined 
as fractions of the light-cone momenta P ± = E±P" of the beam (b) and a target (t) nucleon 
carried by the lepton pair as follows, 



Xi = j^,x 2 = ^L. (6) 



Therefore pj~ t = ap^ = x\ P 6 + , and the coherence length @ averaged over a and Ut reads 



K DY ) 

ir) = - — ^ , (7) 

m N x 2 



where {K DY ) M 2 Xl is 
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In perturbative calculations, one would take m q to be the current quark mass. However, 
this is unsatisfactory in the present case because it leads to large transverse separations 
between the parent and recoil quark when the radiated dilepton takes essentially the entire 
initial quark momentum M 2 (l — a) ~ m q . Indeed, in this case the separation becomes 
r T ~ l/m q , and one can grossly overestimate the contribution of large separations when 
using the approximation cr qq {rT) oc r\. This divergence can be regularized by introducing an 
effective quark mass that suppresses the probability of large separations. Experience with 
DIS |23[] shows that employing m q ps 0.2 GeV gives a good description of DIS data. We use 



the same effective mass in what follows for the DY reaction. Of course, this divergence does 
not create any problem when the dipole cross section saturates at large separations. 

Clearly the coherence length, which controls nuclear shadowing, may not scale in X2 and 
differs from the usually- used approximation l c = l/2x2m^. Moreover, it follows from (|^) 
that as xi — *■ 1, l c — > (since a > x±), i.e. nuclear shadowing vanishes. This fact is at 
variance with the usual expectation that shadowing increases towards the limit X\ = 1, which 
corresponds to the smallest x 2 - Although this contradiction shows that QCD factorization 
breaks down dramatically at large Xi, the result is quite consistent with the observation 
that all physics becomes soft [22| at large X\. 



The averaging in (g) should be weighted by the light-cone wave function of the qll 
fluctuation squared, which is known to diverge (for transversely polarized virtual photons) 
at small transverse q — 11 separation tt, i.e. for large kx |H| (similar to the case of deep 
inelastic scattering (DIS) for the fluctuation 7* — > qq). Such fluctuations with small tt 
do not interact; therefore, to avoid the nonsensical result that (K DY ) = 0, which would 
arise from the fact that the vacuum fluctuations are dominated by large kx - leading to 
a divergent normalization - one also needs to include a qq (rT,s) |L6|, [18[] in the evaluation 
of Eq. (|3|). It is also necessary to include it here since we are interested only in those 
fluctuations that participate in the interaction. The dipole cross section plays the same role 
in Eq. (JTp , namely it prevents the interaction from discriminating between the fluctuation 
\q) and \qll) as Tt — > 0. A distinction between the two Fock components is clearly needed 
in order to liberate the // pair. 

In what follows, we are interested in nuclear shadowing, which in lowest order originates 
from double scattering. For reasons given above, the relevant mean coherence length should 
be weighted with a^ q (ar T , s). Thus, Eq. (|S|) can be written explicitly as, 
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J dx q F } q l (x q ) J d 2 k T m- Uq (x 1 /x qi k T ) K DY {x x jx q , k T ) 

' M 2 ,X! 



k dy ) 9 =™ — — — : : — , (9) 



/ dx q F£(x q ) J d 2 k T ^(xx/xq, k T ) 



xi 



where F q (x q ) is the quark distribution function of the beam hadron, which depends on the 
fraction x q of the hadron light-cone momentum carried by the quark. Hereafter we restrict 
ourselves to the transversely polarized 11 pairs dominating the DY cross section, noting that 
the contribution of longitudinal polarization vanishes in any case as X\ — > 1. The modified 
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light-cone distribution amplitude ^>J lq reads |23 



D ir T -k T 



*B*(«> M = I d 2 r T e"^- a, g (ar T , s) r T ) 



where pi jl^, p] 



^ 9 («,M = ^^X/0 T x^o(rr T ). 



(10) 



1 o ^> 
2 7T 

Here Xi,/ are the spinors of the initial and final quarks, Z q is their charge, K (x) is a modified 
Bessel function, and 

(12) 



[1 - a) M' z + a 2 m 2 



The operator T has the form 

T = i m q a 2 e* ■ (nx a) + a e* ■ (a x V) — «(2 — a) e* • V 



(13) 



where e is the polarization vector of the // (virtual photon), n is the unit vector along the 
projectile momentum of the quark, and V acts on tt- 

To simplify calculations we use the dipole cross section in the form corresponding to 
small tt (a more realistic shape, leveling off at large tt f24], 25| , leads to similar results), 



Vqq(rT,s) = C(s) T T , 



(14) 



where the constant of proportionality C(s) does not enter into the final results. Then, the 
modified distribution amplitude in /^-representation has the form, 



*L(a, k T ) = 2 Z g y/a^C(s) e ■ k T ■ 



i a 2 t 2 



7T T* 



+ k 2 ) 3 



(15) 



and the integration over k^ in (|9]) can be performed analytically, 

2M 2 N DY 



K 



DY 



where 



N~ 



DY 



da FHxi/a) 



3 D DY ' 
l-a) a 2 [l + (l-«) 2 ] 



(16) 



(17) 



D DY = JdaFfcJa) - [1 + ^ 



a 



■t-i 



The coherence length 1® Y calculated with (|7|) and (|T^) - (]18[) is plotted by solid curves as 
function of x\ in Fig. ^|for different dilepton masses M = 4, 5, ... 8 GeV. 

Equation (|16|) should be compared with the analogous factor (K DIS ) = m N x 2 {lc IS ) for 
DIS calculated in Ref. |2"3"[, 

, DI8 \ 2Q 2 N™ 



K 



3 D DIS 



(19) 




0.2 0.4 0.6 0.8 1.0 

/"VB 



Figure 2: The mean coherence length as function of x± and dimuon 
effective mass M = 4, 5, ... 8 GeV (solid curves). Dashed lines show the 
predictions of QCD factorization which relates the cross section of the 
DY reaction and DIS on a nucleus. 



where 

i 

t dis 



a(l - ol)[ol 2 + (1 - a) 21 



N = da ^ (20) 



e 6 



o 

D DIS = i <ln - ■ \ •" . (21) 



a 2 + (l-a) 2 



o 



e 4 



and where 

e 2 = a(l-a)Q 2 + mJ . (22) 

Obviously, (K DIS ) depends only on Q 2 , while (K DY ) depends on M 2 and X\ and, as we 
have noted, vanishes as ii — > 1. 
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We compare in Fig. ^ the coherence length (l c ) = (K)/m N x 2 calculated with (|T6| ) and 
(fl9f) for kinematics corresponding to the E772/E866 experiments. The dashed curves are 
calculated for DIS at the same XBj = x 2 . The discrepancy between the two sets of curves, 
which increases towards x\ — 1, manifests a deviation from factorization. Nevertheless, 
factorization is restored at small x\ and/or at large M. Note that the factorization theorem 
requires only that the soft physics, which is common to DIS and DY, should factor from the 
reaction mechanism at large Q 2 . Since shadowing in the distribution function is controlled 
in the target rest frame formulation by the coherence length and an effective cross section 
{(J e ff discussed in Sect. |] below), factorization requires that these two quantities should also 
be the same for DIS and DY in this limit. 

Note that the x 2 dependence of the mean coherence length is controlled mainly by the 
denominator in (0), while the factor (K DY ) is a rather flat function of x 2 and X\. This is 
demonstrated in Fig. |3] where (1® Y ) is plotted versus x\ for different fixed values of x 2 . 

4 



3- 



2 



E=800GleV 




x 2 =0.08 



0.1 0.3 0.5 0.7 0.9 

Xi 

Figure 3: The mean coherence length (0) as function ofx\ at fixed values 
ofx 2 = 0.02, 0.03, ... 0.08. 



At fixed x 2 the coherence length is nearly constant within the kinematical range of the 
E772/E866 experiments at x\ < 0.8. The x 2 dependence of (K DY ) Xl)X2 at x\ = 0.7, depicted 
in Fig. f|, demonstrates that (K DY ) = 0.2 — 0.3 (see also in p3| ) is quite small compared to 
the usually accepted (K DY ) = 0.5. In the kinematical range of the E772/E866 experiments, 
M > 4GeV, and (K ) ~ 0.3 does not vary much. 
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Figure 4: Factor K DY appearing in the expression for (j^ Y ^ (Eq. as 
function of x 2 at X\ = 0.7. 



3 Energy loss of the projectile quark in nuclear matter 

The coherence among the soft partons of the incident proton is destroyed as a result of 
inelastic interactions of the projectile in the surface of the nucleus. These partons then 
propagate through the nucleus independently, losing energy (to hadronization) and having 
their transverse momentum broadened by interactions. This is predominantly a nonpertur- 
bative process, since the projectile quark contains mainly soft fluctuations. 

Eventually, one of these quarks may develop a hard fluctuation containing the lepton 
pair such that the lifetime of the fluctuation is short compared to the mean spacing between 
bound nucleons. The 11 pair is immediately produced when a hard interaction occurs with 
a bound nucleon within the very short fluctuation time that we have referred to as the 
coherence time. Any change in the longitudinal or transverse momentum of the projectile 
quark caused by the preceding initial state interactions can diminish the energy of the 
projectile quark participating in the DY reaction. This affects the momentum spectrum for 
the produced DY pair, which then serves as a probe for the dynamics of the initial state 
interaction. It corresponds to the first of the two sources of nuclear suppression of Ra/d(%f) 
mentioned in the introduction. 
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3.1 Models for energy loss 
3.1.1 String model 

The first inelastic interaction of the beam in the surface of the nucleus via color exchange 
with a bound nucleon leads to the formation of color strings between the target and beam 
partons. Due to the constant retarding action of the string, the leading projectile quark loses 
energy with a constant rate per unit length, dE/dz = —k s [[], ^BJ, which is invariant relative 
to longitudinal Lorentz boosts. The string tension k s is related to the slope parameter a' R 



of the mesonic Regge trajectories p7| as 



k. = — J-rwlGeV/Zm. (23) 
1 7r a R 

This value imposes a scale for the expected rate of energy loss, which is independent of time. 
Correspondingly, the energy loss increases linearly with length L of the path, AE = n s L. 
The energy lost goes into acceleration of the target quarks and production of new hadrons. 

In the simplest version of the string model, where at most one string can attach to a given 
quark, multiple interactions of the leading quark in the nuclear medium are expected to have 
no influence on its energy loss |28|]. Indeed, no matter what happens to the leading quark 
before the hadronization is completed, it remains a color triplet and is slowed down with the 
same density of energy loss k s by the color-triplet string attached to it: the quark is always 
being retarded with the same k s . Although such a conclusion sounds puzzling, it does not 
mean that multiple interactions in nuclear matter have no effect on hadronization. These 
interactions do make the fragmentation function of the quark softer, i.e. more particles 
are produced in the fragmentation region of the nucleus (dn/drj oc A 1 / 3 ) as the leading 



hadron diminishes its energy H |26j. This indeed may look like the result of energy loss, but 
it happens because each rescattering of the color string in nuclear matter initiates a new 
hadronization process from the very beginning, but with decreased initial energy |28[ . 

The above treatment of the interaction is obviously oversimplified; in fact, multiple 
interactions in nuclear medium do induce extra energy loss. Indeed, the the valence quarks 
are surrounded by parton clouds, only a part of which can be resolved by a soft interaction. 
Multiple interactions in nuclear matter are obviously able to involve more of the partons. 
Correspondingly, the energy loss should grow more steeply than linearly in L. 

This idea is realized in the dual parton model [^J (or quark-gluon string model |3(|) by 
assuming that each multiple interaction in the nuclear medium activates a new sea qq pair, 
leading to the formation of an extra couple of color strings. The probabilities of multiple 



string production are given by the Abramovsky-Gribov-Kancheli (AGK) cutting rules |H[ 



The dual parton model skips over the space-time development of string production and 
decay, which takes long time proportional to the initial energy, jumping directly to the 
final spectrum of produced particles. We, however, are interested in the early stage of 
hadronization when the projectile quark is still propagating through the nucleus and its 
cloud of sea partons is still coherent. Therefore, we should follow the overall energy-loss of 
a nonperturbative constituent quark retarded by more than just one color string. 
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A quark covering the distance z in the nucleus, from the point of the first inelastic 
interaction of the incident hadron, loses an energy per unit length given by 

^ = «.[l + (n(z)>] , (24) 

where (n(z)) = a qN pA z is the mean number of collisions experienced by the quark over the 
distance z. One can approximate the value of the nuclear density by pa ~ 0.16 fm~ 3 , since 
the path between the first inelastic collision and the DY pair production covers essentially 
only the interior of the nucleus. Correspondingly, the energy loss over a distance L acquires 
a correction oc L, 

AE = k s (l + - a qN p A L 2 ) , (25) 



2 

The phenomeno logical quark- nucleon cross section is model dependent. In the additive 
quark model, a qN = a^ N /3 ~ 10 mb. The same value follows from the dual parton model 
|29| , |30|1 , which uses the weight factors <r n for the n-fold scattering appearing in the Glauber 
model as powers of er^ . The new strings formed after each rescattering should be shared 
by three valence quarks. 

In a more realistic model the main part of the hadronic cross section corresponds 
to a soft nonperturbative interaction, which is unable to resolve and free the gluons that are 
located with small transverse separation (tt) = ro = 0.3 fm around the valence quarks [ f25|j . 
This part of the cross section is independent of energy and does not obey quark additivity. 
Clearly, such soft interactions of the quark skeletons of the proton contribute to the linear 
term in (^) . Only a semi-hard interaction is able to resolve the small-size gluon cloud and a 
correspondingly rather small quark-nucleon cross section related to the excitation of quark 
and gluon radiation. This process corresponds to the higher terms in topological expansion 
in the dual parton model, i.e. contributes to the L 2 term in (p5|) . The predicted cross section 
is a qN = (9/4) rl (s/s ) A , where §A = 0.17 and s = 30 GeV 2 . At s = 1600 GeV 2 } one 



gets a qN = 4.2 mb, more than twice as small as the prediction of the additive quark model. 

The mean values (L) and (L 2 ) in a heavy nucleus (with constant density) corrected for 
the mean free path A ~ 2 fm of the incident proton in the nucleus are, 

(L) = \Ra-X 

(L 2 ) = \r 2 a -\r a \ + \ 2 ■ (26) 

For example, for Ra = 6 fm the nonlinear L 2 correction in (|25|) is only 24% and 10% in 
the two models discussed above. It decreases dramatically, of course, for lighter nuclei. We 
neglect this small correction and assume for further applications that energy loss is a linear 
function of L. 

3.1.2 Perturbative QCD 

Although application of perturbative methods to soft processes is not completely legitimate, 
one may hope to get at least the scale of the effect by doing so. Energy loss treated 
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perturbatively originates from gluon bremsstrahlung by a quark propagating through a 
medium. Here again one should distinguish between two sources of energy loss, namely gluon 
radiation originating in the first inelastic interaction, which occurs even in the vacuum, and 
gluon radiation induced by multiple interactions of the quark in the medium. 

It was first demonstrated by Niedermayer |33] that as a consequence of coherence effects, 
(photon) gluon bremsstrahlung in the vacuum carries energy away with a constant rate, 

%max j 

AE(L) = E q J d 2 k T J dxxj^ , (27) 

o T 

where &t and x are the transverse momentum and fraction of the quark light-cone mo- 
mentum of the parent quark momentum carried by the gluon, respectively. The upper 
integration limit x max is fixed by the condition that all gluons contributing to AE are ra- 
diated within the path L. In order to be radiated, the gluon must lose coherence with the 
parent quark, otherwise one cannot disentangle between the quark plus its color field and 
the quark plus the radiated gluon. The distance If over which this loss of coherence occurs 
and radiation formed is specified by the condition that the quark-gluon separation must 
exceed the transverse wave length of the gluon, which leads to If = E q x/k^, where E q is the 
quark energy [j34 |. Note that this value is twice as small as the coherence length defined in 
(@). Applying the condition If < L we get x max = Lk\jE q . 

The spectrum of radiated gluons with k\ smaller than the mean momentum transfer 
squared in the collision was calculated in [[35H . 

(28) 



dn g 3a s 



dx d 2 kx 7r 2 x k\ 



Thus, the density of energy loss reads, 

d E 3a, 



d z 71 



(4) , (29) 



and a quark loses energy in the vacuum with a constant rate proportional to (k^). Note 
that this is not the induced energy loss in a medium. 



The proper value of a s ~ 0.4 for gluon radiation was calculated in ||32| . The transverse 
momentum squared of radiated gluons is rather large kx ~ 1/to = 650 Me V [25] , |32|. 
With these values we evaluate the energy loss rate (p9|) as dE/dz ~ 0.8 GeV/ fm, which is 
amazingly close to the string tension. Note that the radiative energy loss is not supposed to 
be an alternative to what is given by the string model, but it is rather a different source that 
should be added to it. This is analogous to the additive contributions to the total hadronic 



cross section arising from nonperturbative strings and gluon radiation [32 

Quite a different situation arises if the quark originates from a hard reaction as, say, 
in DIS or a high-p^ scattering event, where the rate of energy loss from gluon radiation is 
expected to be tremendous 0. 

Another source of energy loss is induced gluon radiation due to multiple quark interac- 
tions in the medium. This has received much attention recently P?L Kfl], since it may serve 
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as a probe for the production of dense matter in relativistic heavy ion collisions [36|. The 
induced gluon radiation was calculated perturbatively |36|], and it was concluded that the 
energy loss grows like AE oc L 2 . This is a direct consequence of the relation ( |29|) between 
the rate of energy loss and the transverse momentum squared of gluons, which follows that 
of the quarks and is known to rise oc L. The relation between induced energy loss and 
nuclear broadening of the mean transverse momentum squared of the quark was established 
in [ |36f| , 

AE=^A(&)(L), (30) 

where the broadening of (p^) was measured J57I] to be rather small A(p^) « 0.1 GeV 2 even 
for tungsten. Therefore the coefficient of (L) in (|30|) is about 0.075 GeV/fm, an order of 
magnitude smaller than the string tension k s . 

The induced energy loss of a nonperturbative quark, both the string and radiative parts, 
turn out to be a rather small fraction, less than ~ 10 — 15 %, of the total energy loss even 
for heavy nuclei. Although this correction has a different L-dependence, one can effectively 
absorb it into the main linear term. 

The dominant constant parts of the energy loss rate dE/dz (so-called vacuum energy 
loss), related to the string tension (Sect. 3.1.2) and to gluon radiation, have different origins 
and, as we have said, should be added. The resulting rate of energy loss is thus expected to 
be 

-4^ ^IGeVjfm . (31) 
dz 



3.2 The quark path in the nucleus 

The DY reaction on a nuclear target is usually mediated by the debris resulting from an 
inelastic collision between the incident hadron and a bound nucleon on the front surface of 
a nucleus. This debris, once produced, propagates through the nucleus and produces the 
observed lepton pair when it strikes a bound nucleon, as illustrated in Fig. [5] (top). One 
would not expect to observe any difference between this and the DY reaction on a free 
nucleon for z 2 — Z\ — > 0. The reason is that the primordial momentum distribution of the 
projectile partons cannot be affected by soft inelastic interactions occuring over short times. 
However, when the distance z 2 — Z\ is finite, the soft projectile partons may lose energy 
through hadronization and thus reach the point z 2 with diminished energy. However, at 
fixed Xi the DY pair has to be produced with the same measured momentum, i.e. with 
an increased fraction of the initial momentum. As a result, the cross section for a DY pair 
produced with a given longitudinal momentum turns out to be smaller on a nucleus than it 
is on a free nucleon target. 

It is usually incorrectly assumed that the quark propagates from the surface of the 
nucleus to the point where the DY pair is produced, which would mean that the mean quark 
path in the nucleus would be (L) ~ 3i2^/4. According to Fig. [|, this should be shortened 
by at least the mean free path of a proton in a nucleus, 2 fm. This would substantially 
reduce (L), by a factor of two or more, so that the mean path between the point of DY pair 
production and the first inelastic interaction is actually shorter than the maximum possible 
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Figure 5: The space-time pattern for DY pair production off a nucleus. 
The upper example illustrates when the beam hadron experiences 

a soft inelastic interaction prior the hard interaction in which the 11 pair 
is produced. The case where the DY pair is produced in the first inelastic 
interaction is illustrated at the bottom. 



distance to the edge of the nucleus shown in Fig. [|. Additionally, there is some probability 
(dominant for light and medium-heavy nuclei) that the incident hadron has no interactions 
prior to DY pair production at point z 2 . In order to find the mean path length (L) of the 
projectile quark in nuclear matter we should average [z 2 — z\) in accordance with Fig. |, as 
done in [Dfl, 



(L) 



0) 



A 



d 2 b / dz 2 p A (b,z 2 ) / dzip A {b,z l )(z 2 - zi) 



x exp 



7 -i 



hN 



dz p A {b,z) 



(32) 



The exponential factor requires that there is no inelastic interaction of the beam hadron 
prior point Z\. The probability of no inelastic interaction of the beam hadron in the nucleus 
prior the DY reaction (see the bottom part of Fig. |^) Wo reads, 

MA 



W 



A a hN 



1 _ TQ>) 



A 'in^ 



(33) 



We calculated (L) for a pA collision taking a. 
parameterization for the nuclear density [^S 



pN 
in 



Pa{t) = p A 



where 



Pa 



1 + exp 
3A 



4ttR a 



1 + 



2 2 \ ~ 1 



n A 



30 mb and using a Woods-Saxon 

(34) 

(35) 
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and the nuclear radius Ra and the edge thickness c are fixed for each nucleus by a fit to 
data on electron- nucleus scattering . 

As expected, (L) is much shorter than 3Ra/4. Both of these results are plotted versus 
A 1 / 3 in Fig. || The corresponding probability distribution in L is given by the expression 



e-i — i 




12 3 4 5 6 



A 



1/3 



Figure 6: The mean length of the path of a quark between the points of 
the first inelastic interaction of the beam proton and DY pair production 
calculated with Eq. (|3^) as a function of A 1 / 3 (solid curve). The dashed 
curve shows the usual expectation (L) = 372^/4. 



W(L) = W S{L) + W^L) 



(36) 



where 



Wi(L) 



r hN 



J d 2 b J dz 2 pA{b,z 2 ) I dz 1 p A (b,zi)8(z 2 - z x - L) 

-co 



x exp 



-a™ J dzp A (b,z) 



(37) 



which is normalized as f£° dL Wx(L) = 1 — Wq. 



Although they are not as large, the effects of energy loss affect the DY cross section on 
a deuteron target as well. In this case the constant is given by the standard Glauber 



formula 39 



NN 
1 _ tot 
16 7T 



J dq 2 F D (4q 2 )e 



(38) 



where B NN ~ 1 ° ^"i/^ 2 ^--i x?_(„i 

i 

fitted in EDf to electron scattering data as 



el ~ YlGeV 2 , and F^q 2 ) is the deuteron charge formfactor parameterized and 



F D (q 2 



Cl e~ biq2 + c 2 e- b2 " 2 



(39) 



Here d = 0.55, c 2 = 0.45, &i = 19.66 GeV" 2 , and b 2 = 4.67 GeV- 2 . This leads to = 
0.96. Correspondingly, the distribution function W\{L) for a deuteron is given by the wave 
function squared related to the formfactor (|39|), 

L 2 f ci 



+ 



c 2 

(4 6 2 + J B^) 3 / 2 



(4 6 1 + 5 e ^) 3 / 2 
/ L 2 



exp 



4(46i + ^ 



cxp 



4 (4 ft, + 2^ 



(40) 



3.3 Nuclear suppression caused by energy loss 

Since the DY cross section decreases steeply as (1 — Xi) n for xi — > 1, any loss of energy by 
the quark by the time it reaches the point z 2 will result in a suppression of the DY cross 
section at large X\. This effect was first suggested and estimated in 0. In the more detailed 
approach taken here, we want to focus on the propagation and energy loss of the projectile 
quark that takes part in the DY reaction. In this case, the ratio Ra/n of the DY cross 
sections of a nucleus to a nucleon can be represented as 

oo 1 

fdLW 1 (L) j dx q F q h (x q )da q D N Y (xl)/dxl 

Ra/nM = W + 1 (Xq)m " 1 . (41) 

/ dXqF^Xq) d u 9 oy{xI) / 'd x\ 

XI 

Here F^{x q ) is the quark distribution function in the incident hadron; x q and x\ = X\jx q 
are the fraction of the light-cone momentum of the incoming hadron h carried by the quark 
and the fraction of the quark momentum carried by the 11 pair, respectively. The lower 
integration limit is given by (x q ) m i n = x\ +AE/Eh, and x\ = x\j{x q — AE/Eh). According 
to the results of Sect. 3.1 we assume that the rate of energy loss dE/dz = —n is constant 
i.e. AE = kL. Note that as the cross section d<7 q DY {xi) / dx\ in Eq. (|4l|) corresponds to an 
incident constituent quark, while that in (|1|) corresponds to a current (perturbative) quark. 
To emphasize that these are different quantities, we use variable x\ in (|4l|) instead of a as 
in®. 

The energy loss AE varies due to fluctuations in the string tension, the number and 
energy of radiated gluons, and due the to the dependence of the distance z 2 — Z\ shown in 
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Fig. |5| on the distribution W(L). We neglect the variation of k, assuming k = const, and 
integrate the numerator of (f|l|) over AE weighted by the distribution W(AE), with W(L) 
given in (|36|). 

The first term on the r.h.s. of Eq. ([41]) corresponds to the first term in (|36|) , which 
represents the probability of no interaction preceding the DY reaction (the bottom part of 
of Fig. |5|). The second term corresponds to one or more projectile interactions as illustrated 
in the upper part of Fig. [| 



3.3.1 Distribution of the valence quarks in the proton 



As mentioned above, in the limit of short coherence length <C Ra, the radiated photon 
arises from a hard DY fluctuation (one lasting for a short time) in the projectile quark, 
whereas the energy-loss mechanism is dominated by the softer fluctuations. Therefore, to 
calculate ( fill) we need to know the distribution function FHx q ) and, correspondingly, the 



DY cross section d o~1^(xl) / d x\ for soft quarks. 

One should rely on a phenomenology of soft hadronic interactions to find the soft quark 

The proper approach is found in the dual parton model 



distribution function F„ 



, X n 



29J (quark-gluon string model |30[|), which describes inelastic hadronic interactions via 
convolution of the quark distribution functions with the fragmentation functions of quarks 
to hadrons. The success of this phenomenology in describing the available experimental 
data on soft hadronic interactions justifies the quark distribution functions used above. 



The proton quark distribution function in 
up and down quarks, 



Fit?* 
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Fu( X q 



If) includes contributions from the valence 



(42) 



where according to [29, 30 



KM = n 



I H I — Xq 

5 8nBPP y q 



(43) 



F p d i 



X n 
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6 (1 



x, 



\2.5 



1 



8 rrP N 
° a tot 

7 8-jtBpp 



(44) 



Here is a normalization factor inessential for (4"ID; the factors 6/5 and 8/7 result from the 
normalization condition as well. The second terms in the square brackets in ( fi3|) and (|4~i|) 
are the first unitarity corrections (we neglect the small higher-order terms). They contain 
the total pp cross section af^ = 40 mb, the slope of elastic pp scattering, B pp = 12GeV~ 2 , 
and an extra factor of 2 dictated by the AGK cutting rules [31]. The behavior of these 
distribution functions — > and cLS Xq ^ 1 is dictated by Regge phenomenology and 
is related to the intercept of the leading meson trajectories ckr(O) = 0.5 and the nucleon 
trajectory a N (0) = -0.5 p9|, BD|. 
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3.3.2 Quark-nucleon DY cross section 



In order to evaluate R e jj s ^{xi) in Eq. (41) we first of all have to fix the x\ dependence of the 
quark-nucleon DY cross section. As was emphasized above, in the limit of short coherence 
length the quark that experiences initial state interactions is a constituent quark requiring 
the distribution function for soft quarks. Therefore, the DY quark-nucleon cross section 
daf^/dxl should be treated correspondingly. Its shape cannot be predicted reliably since it 
depends in an essential way on the model for constituent quarks, which lies in the domain 
of nonperturbative effects. Instead, we parameterize the cross section in the form, 



d 2 a q D N Y 
dM 2 dx1 



K{M 2 )(l-x(f , (45) 



where K(M 2 ) and m are fitted parameters. We use data ^1| for the DY cross section in 
p + 2 H collisions and fit it by the expression 



drr NN \ 
aG DY _ / > jphi 



dxi 




oc dx q F q h (x q ) 1-- , (46) 



with m = 0.362 ± 0.027 taken to be independent of M 2 . An example of the fit to the E772 
data on deuterium target is depicted in Fig. [7| for the dimuon mass interval 7 — 8 GeV. 
effects, 



4 Nuclear shadowing: DY versus DIS 

In the target rest frame, it is clear that the longitudinal momentum transfer to the nucleus 
in the DY reaction becomes small at high energies, in spite of the large effective mass M of 
the dilepton. Therefore, different nucleons start interfering (destructively) in the production 
of the DY lepton pair, a phenomenon usually called shadowing. 

Naively, one would not expect significant shadowing for hard reactions such as DIS, DY 
processes, heavy flavor production, etc., because these processes have tiny cross sections and 
because bound nucleons clearly do not shadow each other. However, substantial shadowing 
does exist for all of these, particularly the DY process, at high energies. The explanation 
is simple. First of all, shadowing arises from the soft components present in the hard 
reaction |16 |. These soft components are small because the probability to develop a soft 



fluctuation in a hard reaction is small. However, this small size is compensated by a large 
interaction cross section. Once a soft fluctuation is created, the reaction is then driven by 
the subsequent hadronic interactions with the nucleus; these interactions are distinct from 
those producing the fluctuation to begin with and may be significant when the lifetime of 
the fluctuation (coherence time) is long. 

The light-cone dipole representation of this reaction is especially suitable for the cal- 
culation of nuclear shadowing. The DY cross section in this representation has the form 
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Figure 7: The cross section of the DY reaction on deuterium in arbitrary 
units as function of X\. The data are from the E772 experiment for 
dimuon mass interval 7 — 8 GeV. The curve is the result of the fit oc 

(l_ Xl )0.362_ 



[16, 17, 18 
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ao NN 

dM 2 dxi 
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da 



a- 



d r T ^n q (a,r T ) a qq {ar T ) 



(47) 



where ^i lq (a,r T ) is defined in (|Tl"l). 

If l c 3> Ra-i the gZZ projectile fluctuations are frozen by Lorentz time dilation during 
their propagation through the nucleus. Therefore, these fluctuations are eigenstates of the 
interaction and should have no inelastic corrections. One can simply eikonalize the dipole 
cross section in (H7J) in the case of a nuclear target N2TL IIBL |18| , 



(Jq q {ar T ) 



2 d% 1 



exp 



^aq q (ar T )T(b) 



(48) 



where the averaging over a and tt is performed in accordance with fl47|) . In the case of weak 
shadowing, one can expand the exponent, dropping the higher-order shadowing corrections 



22 



to obtain 

R A/N = "Ay = 1 " 7 *«// ( T ) + <>{(o>)) » ( 49 ) 



'AT 

where 

(o-| g (ar T ) 

(Teff(Xl, X 2 ,S) = -) ? r{- , (50) 

and where 

(T) = \\ d%T\b) (51) 
is the mean value of the nuclear thickness function 

CO 

T{b) = f dz PA (b,z) , (52) 



with pa{t) the nuclear density distribution. 

Note that naively one may expect a very small transverse size ~ 1/M 2 for the \qll) 
fluctuations that define the dipole cross section in P8| ) and cr e ff- Although these fluctuations 
appear with a vanishing probability at large M, soft fluctuations with a — > 1 can have a 
large transverse separation [[16], [see (|TT| ) - (|i2|) 1 



as follows from (O). For this reason, these soft fluctuations have a large interaction cross 



cross section and can make a sizable contribution to cr e ff, in contrast to the harder ones [[16 



The transition region between the no-shadowing (l c <C 2fm, or X2 ~ 0.1) and asymptotic 
regimes (Z c 3> (L), or l/a;2 3> 3mjviiU), where Eq is applicable, is most complicated 
and needs sophisticated calculations based on the path- integral approach [[43], |23|. Never- 
theless, if shadowing corrections are small (which is the case for the E772 and E866 data, 
where we face only the onset of shadowing), Eq. (pE9|) can be easily interpolated [[42], 



Rf/N « 1 - 7 ^ff < T > F a(1c) , (54) 



4 

where F\{q^) is the longitudinal formfactor of the nucleus 



dze iqcZ p A (b, z) 



2 



(55) 



Strictly speaking, this formfactor must also be involved in the averaging procedure for a e ff 
since q c = l/l c depends on a and kx- In order to speed up the fitting procedure (see below) 
we replace l c by its mean value (l c ). Comparison with exact calculations for DIS done in 



Ref. [E3T] demonstrates that this approximation is sufficiently accurate. 
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Beryllium is too light for the Woods-Saxon parametrization ( 541) of the nuclear density 
with which we describe heavier nuclei. Instead, we use for beryllium the harmonic oscillator 
density |3j 



^2 \ / 



pBeir) =p° Be \1 + a — ] exp ( -— ] , (56) 

where 



a 2 I \ a 2 



PBe = , 7^ . (57) 



i.Tra 2 ) 3 / 2 (l + fa 



a = 1.77 fm, and a = 0.631. 

The expression Eq. fl5"5|) is designed for medium and heavy nuclei and cannot be applied 
to the deuteron target involved in our analysis. In this case, one should use a different 
expression |39j [compare with (|38|)1, 



oo 



Rd/n = 1-J^J dq 2 T F D (4q 2 ) , (58) 
o 

where q 2 = q\+q 2 c , with q? the transverse momentum. Using the deuteron charge formfactor 
F]j(q 2 ) in the of form Eq. (|39"D, we arrive at shadowing for a deuteron, 

The effective cross section (a 2 ) / (a) can be calculated by averaging in accordance with 

^ = £ daFP( Xl /a) [1 + (1 - a) 2 ] r 2 /q 2 J d 2 r T K 2 (rr T ) <rj,(arr,s 2 ) 
e// 1 ' 2 ' daFq(xifa) [1 + (1 - a) 2 ] r 2 /a 2 J d 2 r T K 2 (rr T ) aq q (ar T ,x 2 ) 

where Ki(y) is the modified Bessel function and r is defined in Eq. (|T2"D. It is reasonable to 
use here the phenomenological x-dependent dipole cross section |24|] which describes very 
well data from HERA for the proton structure function at high Q 2 , 



Vqq(P,x) = O"0 



( P2 

1 - exp j— 



(61) 



where <To = 23.03 mb and po(x) = 0.4 (x/xo) OAM fm, xo = 3.04 ■ 10~ 4 . In the case of the DY 
reaction, x = x 2 . 

The results of the calculations for a e ff(xx,x 2 , s) at s — 1600 GeV 2 are depicted in Fig. || 
as a function of x\ for different fixed values of x 2 . Note that the effective absorption cross 
section a e ff substantially increases with x\. This is a manifestation of factorization breaking, 
since suppression {i.e. shadowing) of the nuclear structure function cannot increase with 
rising M 2 at fixed x 2 . In the light-cone approach this behavior is easily interpreted: the 
mean transverse separation squared that enters into (p0"D, a 2 (rfi) ~ 1/r 2 « 1/[M 2 (1 — a)], 
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Figure 8: The effective cross section Eq. ^) as a function of x\ for 
different fixed values of x 2 = 0.02, 0.03, 0.05 and 0.07. 



increases as x\ — > 1 because a > X\. Nevertheless, Fig. || shows that the x\ dependence of 
a e ff is quite weak for xi < 0.8. 

Although we cannot test our calculations of shadowing by comparing them to nuclear 
suppression in DY reactions, which are affected by energy loss effects, we may check them 



against nuclear shadowing in DIS. In this case a e ff reads |23 



Dis ( q2 ) = Jo da [a 2 + (1 - a) 2 } e 2 J d 2 r T K 2 {er T ) a 2 q {r T , x 2 ) 
r eff [x 2 , Ui ) £ ^ [a2 + (1 _ a)2] £2 J K ^ erT) ^ 



where e is introduced in ([22]). 

The calculated a^j I ^{x 2) Q 2 ) turns out to be rather different from a^f(x 2 , Q 2 ) calculated 
with Eq. (|60|). Their ratio is plotted in Fig. |9] as a function of X\ for different dilepton masses. 
The fact that the ratio is not unity is a deviation from factorization (except near x\ = 1, 
where the factorization theorem does not apply since, as we have noted, all physics becomes 
soft here |p2] ). However, the ratio may be shown to approach unity logarithmically for large 
M 2 (and x\ < 1), as required by the factorization theorem and as suggested by the results 
in the figure. 

The effective absorption cross section controlling nuclear shadowing is related to forward 
diffractive dissociation ||44j| . As a simple test, we can compare our results for with data 
for diffractive dissociation of highly virtual photons which were measured at HERA to be 
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Figure 9: Ratio a^/a^/f as a function ofx\ for M = 4 (solidcurve) , 20 
(dashed) and 100 GeV (dotted). Deviation from unity demonstrates 
breakdown of QCD factorization (see text). 



about 10% of the total DIS cross section 



CT DIS CT 



,ff(x 1 ,X 2 , s) 



n DIS i R _ p>DIS 

"tot LG7rn dd 



0.1 , (63) 



where the subscript dd means diffractive dissociation, and B^J 8 ~ 5GeV~ 2 . 

A more rigorous test would involve a direct comparison of the calculated nuclear shad- 
owing for DIS with data. Shadowing can be calculated using Eqs. (|62|) and (|54|) , where 
the argument of the formf actor, q c = xm N (K DIS ), is given by Eq. ( |T9"D (see Fig. g). The 
results are plotted in Fig. [10] as the dashed curve and compared with data from the NMC 
experiment. Since the value of Q 2 in the data varies with x we incorporate this correlation 
in our calculations; we use the same set of x and Q 2 as in the NMC data. Shadowing for 
the deuteron is also taken into account. At large x ~ 0.1, the data display about 2% anti- 
shadowing, the origin of which is poorly understood, although it is probably related to the 
same nuclear medium effects that cause the EMC effect. As this dynamics is not included 
in our calculations, it is not surprising that we underestimate the data in this region. 

At smaller x (and smaller Q 2 ), as expected, the calculated shadowing keeps decreasing 
and does not seem to saturate. This is due to the rise of a e ff(x, Q 2 ) at small x and Q 2 . 

At even smaller x, gluon shadowing also becomes important. The onset of gluon shadow- 
ing occurs only at x < 10~ 2 because it is related to fluctuations containing extra gluons that 
are heavier than qq fluctuations. Correspondingly, the coherence length for gluon shadowing 
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Figure 10: Normalized ratio of carbon to deuterium structure functions. 
The data are from the NMC experiment p5|/ . The dashed curve is calcu- 
lated using Eqs. (]5^), (^) for each data point at the same values 
of x and Q 2 . The solid curve includes also the effect of gluon shadowing 



as it is calculated in |2^/. 



is shorter than it is for quarks p3| . Nuclear shadowing for gluons at small Q 2 relevant for 
the NMC data and at low x is predicted in |25| (see Fig. 6 of that paper). We multiply 
F£ I by the gluon shadowing factor as calculated in [^5] . The result shown by the solid 
curve is in good agreement with the data and confirms the reliability of our description of 
shadowing, which contains no free parameters. 



5 Energy loss versus shadowing 

By including shadowing explicitly in the analysis of the DY data, we are able to utilize 
the entire data set of the E772/E866 experiments for the ratio R exp (x,M 2 ). This must 
be done recognizing that shadowing and energy loss are complementary sources of nuclear 
suppression. If the lifetime of a fluctuation containing the lepton pair substantially exceeds 
the size of the nucleus, we assume that the energy loss of the hadronic part of the fluctuation 
does not affect the spectrum of dileptons, which are created as a fluctuation long in advance 
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of the nucleu^]. Thus, in the two limiting cases of very long or very short coherence time 
one finds either shadowing or energy loss present, correspondingly, not both together. 

The DY amplitude in the long l c limit differs from that in the short l c limit by addi- 
tional terms that account for interference between the DY amplitudes involving different 
numbers of nucleons. Within the approximation that the double scattering correction dom- 
inates, which is quite accurate if shadowing is small (appropriate for the kinematics of the 
E772/E866 experiments), the interference term always appears multiplied by the nuclear 
formfactor squared, as in Eq. (|54|). The transition from the short l c regime, described by 
Eq. (fill) , to the regime of l c 3> Ra, described by Eq. ([54]), is likewise controlled by the 
nuclear formfactor. It is easy to write the interpolating expression, namely 



R A/N (x 1 ,M 2 )= l-F 2 A {q c ) R E A % s { Xl )-l + Rf%(x 1 ,M 2 ), (64) 



L A/N 



->shad( ^ /\/f2\ 



h A/N\ 



where F\(q c ) is given by (|55|) and Ra/°n s an d Ra/n are gi yen by (|4lD and (|54]) respectively. 
Obviously, (Q) correctly reproduces the short and long coherence length limits and mixes 
these effects when both contribute. A nice feature of Eq. (|64|) is that the effect of energy 
loss is weakened with rising l c and eventually vanishes at l c R A . 



5.1 Results 

With Eq. (|6~4]), we have adjusted k to fit the entire set of ratios C/D, Ca/D, Fe/D, W/D, 
Fe/Be and W/ Be from the E772 [T(J and E866 |13| experiments double-binned in x\ > 0.3 



and M > 4 GeV. Since the normalization of the data is subject to systematic uncertainties, 
we included them in the fit by introducing for each experiment an overall normalization 
N = 1 ± AN for the theoretical values of the nuclear ratios. We allowed these to contribute 
to x 2 by treating them as additional experimental points and assuming a normal distribution 



for them. A few examples of the fit results are depicted in Figs. [LI], [12], and [13] (more bins 
in M and a greater variety of nuclei were involved in the fitting). In addition to the solid 
curves, which show the full result of the fit, we also show the contribution of shadowing 
alone in the dashed curves. Shadowing was calculated for the mean value of mass calculated 



for each interval as J (M 2 ). 

We found the following value for the rate of energy loss, 

- ^ = 2.73 ± 0.37 GeV/fm , (65) 
d z 

with x 2 /d.o.f. = 0.9, and normalization factors N(E772) = 1.014 ± 0.005 and iV(£866) = 
0.996 ±0.006. 



lr This is different for nuclear broadening of the transverse momentum of a DY pair. Since a nuclear target 
supplies a stronger kick to the quark part of the fluctuation compared to a nucleon target, the nucleus is 
able to break up the coherence of harder fluctuations, i.e. fluctuations with larger transverse momentum of 
the lepton pair frL~8j. 
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Figure 11: Examples of ratios of the DY cross sections on carbon to deu- 
terium as functions of x\ for various intervals of M. Dashed curves cor- 
respond to net shadowing contribution, solid curve show the full effect 
including shadowing and energy loss. Data from Refs. flK)/,p^/. 



A nice feature of the fit is the consistency of the values of dE/dz resulting from inde- 
pendent fits to the E722 |§ and E866 [|7]] data, 



dE 

dz 
dE 

dz 



£772 



£866 



2.32 ± 0.52 GeV/fm 
3.14 ±0.53 GeV/fm 



(66) 
(67) 



with x 2 /d.o.f. = 0.85 and x 2 /d.o.f. = 1.02 for the E772 and E866 data respectively, and 
normalization factors N(E772) = 1.01 ± 0.006 and N(E86Q) = 1.00 ± 0.007. 

The rate of energy loss given in Eq. (p5|) is larger than the value 1 GeV/ fm suggested 
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Figure 12: The same as Fig. [II], except that the ratio of tungsten to 
deuterium is shown. 



by the string model. However the string tension, n s = 2 7ca' R ~ lGeV/fm is a static 
quantity, related to the slope parameter a! R describing the orbital excitation spectrum of 
hadrons. Evidently the energy loss by a quark includes an additional piece, such as the 
dynamical contribution arising from gluon radiation ||, 26 1. The magnitude of the latter 



was determined in section |3.1.2j to be similar in size to the energy loss arising from the 
formation of strings. Adding this contribution, the result ( |6~5l) agrees with the expectation 
Eq. (0). 

The normalization factors found for the E772 and E866 data are well within the quoted 
systematic errors. These values of norms explain why the dashed curves representing net 
shadowing in Figs. [U], [12| do not match unity for short coherence length. At the same time 
the dashed curves slightly rise above 1 at large x\ and M, demonstrating antishadowing. 
This is a result of the delayed onset of shadowing toward small X2 for heavy nuclei (W) 
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Figure 13: The same as Fig. [II], except that the ratio of tungsten to 
beryllium is shown, with data from Refs. |J3|/, 



compared to light ones (D) as dictated by their formfactors. 

One can see from Figs. [TT] - [L3| that the effects of energy loss and shadowing display 
quite different behaviors as a function of A, x%, and M. For carbon nearly all nuclear 
suppression, to the extent that it exists, comes from shadowing. This is because the mean 
path length available for energy loss in nuclear matter vanishes for nuclei as light as carbon, 
as one can see from Fig. |6|. At the same time, for tungsten, energy loss makes a substantial 
contribution to nuclear suppression, while shadowing vanishes for large masses because the 
nuclear formfactor falls steeply at large x-x- This difference in the A and M dependence 
helps to disentangle the two effects rather effectively. 

Another specific signature of energy loss is a suppression which does not vanish even at 
small x\. Indeed, this effect of the order of 4% seems to be supported by data for tungsten 
in Figs. O, O, while data for carbon, Fig. O show no deviation from 1. This feature also 
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helps to single out the effect of energy loss since the data have especially high accuracy in 
this region of small x\. 



Notice that the ratios of tungsten to beryllium shown in Fig. [13] reveal an interesting 
feature; namely, at large dimuon masses this ratio rises with x\ reaching values above one, 
reminiscent of antishadowing. Instead of this, one might have expected a stronger shadowing 
for a heavier nucleus. Of course Eq. ( [54]) cannot lead to any antishadowing. However, the 
onset of shadowing with decreasing x 2 must be delayed for heavier nuclei whose formfactors 
decrease more steeply with q c . We illustrate this feature in Figs. |14] and [15|, corresponding 
to the dimuon masses M = 4.5 and 7.5 GeV, respectively. We compare nuclear shadowing 
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Figure 14: Shadowing in DY reaction on carbon, iron and tungsten as 
function of x 2 at M = 4.5 GeV. Nuclear shadowing disappears at large 
and small x 2 because the coherence length, Eqs. (0) and vanishes in 
these limits. 

Figure 15: The same as in Fig. [1^, but at M = 7.5 GeV. 



calculated with Eq. (p4[) for carbon, iron and tungsten as function of x 2 . Indeed, one can 
see in Fig. |TJ] that the onset of shadowing for heavier nuclei occurs at smaller x 2 (compare 
with Figs. 5 and 6 in f25|). At the larger dimuon mass of M = 7.5 GeV in Fig. [15|, the 
behavior is seen to be even more complicated. In this case l c is so short that it suppresses 
the formfactor of tungsten at all x 2 . In any case, these tiny variations of shadowing have 
no influence on our final results for energy loss. 

Note that shadowing vanishes towards the kinematical limit x\ = 1, i.e. minimal values 
of x 2 . This is the result of the coherence length vanishing in this limit (see Fig. |||). This 
property is irrelevant for our analysis since there is no data in this region. 
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5.2 Stability of the solution 



We have performed several tests of the stability of our results given in (p5|) - (p7|). First of 
all, as already mentioned, the analysis could be affected by any overlooked physical effects 
that are related to nuclear structure and that exist at large where no shadowing is 
expected. To test this, we varied the relative number of shadowed and unshadowed events 
by imposing an upper cut-off x™ ax on the values of x 2 allowed in the data set. As this cut-off 
is lowered, more points are affected by shadowing. On the other hand, decreasing the cut-off 
x max diminishes the influence of any missed physical effects, as mentioned above. We have 
plotted the values of the energy-loss rate n resulting from separate fits to the E772 and E866 
data as a function of x™ ax in Fig. ITB. The results of the fits to both sets of data appear 
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Figure 16: The rate of energy loss —dE/dz = function of the 

upper cut-off x 2 < x™ ax . Closed and open circles correspond to fits to 
the E772 and E866 data respectively. 



reasonably stable, confirming the correctness of our calculations for nuclear shadowing of 
the DY cross section. 

Secondly, we tested the sensitivity of our results to the magnitude of the shadowing in 
the analysis. We replaced 1 — R^jp C s h a d({ — R s a/d)i where the factor C s h a d was varied. 
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Eliminating shadowing entirely by fixing C s h a d = 0, we found —dE/dz = 2.34±0.37GeV/ fm 
with x 2 — 240 from the common fit of the E772 and E866 data. Next, making another 
extreme assumption, we doubled the amount of shadowing by fixing C s h a d = 2 and found 
—dE/dz = 3.00 ± 0.37 GeV/ fm with x 2 — 259. Whereas one might be tempted to conclude 
from this exercise that the treatment of shadowing is unimportant the extraction of energy 



loss, one clearly sees from from Fig. [16] that this is not true: the energy- loss value determined 
by the data acquires important contributions from the shadowing region. We also checked 
the optimal value of C s h a d by treating it as a free parameter and found C s h a d = 0.88 ± 0.18 
with x 2 = 218, which confirms correctness of our calculations for shadowing. Thus, we 
conclude that by varying the amount of shadowing over a wide range, the rate of energy 
loss remains unchanged within error bars, while x 2 becomes substantially worse. 

Next, we examined the accuracy of the double-scattering approximation in Eq. ([54]). 
We had assumed that if the shadowing effect is weak, we could neglect the triple and 
higher order rescattering terms. The exact calculation of these corrections in [^, [L9| is 
too cumbersome to be used in the fitting procedure. Instead, one can include the multiple 
interaction corrections in the same way as in Ref. [11], i.e. assuming that rescatterings 
happen with the same cross section a e ff that governs the shadowing correction in (|54]). 



<-'2 



Then, the formfactor Eq. ( [5"5] ) should be modified as, 

F a(Qc) = -j=r j * d 2 b J dz 2 p A (b,z 2 ) J dz x p A (b, z{) exp iq c {z 2 - - ^ J dz p A (b, z) 

— oo — oo L z\ 

(68) 

We repeated the fit replacing F\(q c ) =>- F\(q c ) in Eqs. ( pjl ) and (0) and found no visible 
deviation from the result in Eq. (|65"D. 

It was also of interest to know whether our results would change if, instead of using 
Eq. ( |64"D to combine the effects of shadowing and energy loss, we make the maximal-strength 
mixture of these two effects as assumed in |0|; namely, Ra/n = Ra 1 /°n s ^a/n- Fitting all data 



when we do this, we find nearly the same value of the energy loss rate that we previously 
obtained, -dE/dz = 2.18 ± 0.31 GeV/ fm. 

One can also extract from the results of these modifications a scale for the model de- 
pendence of the systematics in the rate of energy loss. We evaluate such a systematic error 
at about ±0.5 GeV/ fm. Note that although all these variations of the fitting procedure do 
not affect the results very much, they do substantially increase the value of x 2 compared to 
the fit done in Sect. [O. 

The pleasing stability of the results confirms the conclusion of Sect. [TTJthat there should 
be no confusion between the effects of shadowing and energy loss (provided that the data 
are sufficiently exact and binned both in x\ and M) since these two effects behave very 
differently as functions of A, x\ and M. 

Finally, we would like to call attention to the fact that the data we are analyzing for k 
cover the kinematical region of the target x extending all the way out to large values of x 2 . 
This is a concern, since for x 2 > 0.1 there exist EMC and antishadowing effects in nuclear 
quark distribution function that we have not yet accounted for in the theory. We suppress 
the contribution of large x 2 selecting data with x\ > 0.3. Besides, we argue next that these 
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effects make a negligible contribution to the data we are analyzing. 

To show this, we first remark that the EMC/antishadowing effects appear in the quark 
distribution function of the nuclear target, corresponding to antiquarks in the proton pro- 
jectile. We expect that the EMC/antishadowing corrections, if appreciable, can affect the 
analysis only at small x\, since antiquarks in the beam contribute only for small x\. These 
corrections enter the theory as expressed in Eq. fl64|) in both R^0(xi,M 2 ) (Eq. fl54|) ) and 
R E A % s (x x ) (Eq. ©). 

Consider first the shadowing contributions. Shadowing is appreciable only for very small 
x 2 < 0.03. Taking s = 1600 GeV 2 , these small values of x 2 correspond to values of xi = 
M 2 (sx 2 )~ 1 > 0.3 at the smallest values of M in our data set. For x\ > 0.3, the antiquark 
distribution function is quite small, and in any case cr e ff, Eq. (0), is not very sensitive to 
how we average. So, we can neglect the antiquarks in the shadowing region. 

Consider next the energy-loss contribution. For E772, all points in the data set cor- 
respond to X\ > 0.3. This means that the data points are uncontaminated by EMC and 
antishadowing effects (i.e., x 2 < 0.1) for values of M < 6.9 GeV. However, for larger values 
of M in E772, and for all values of M in E866, some region of X\ will be affected. We have 
made quantitative estimates of the size of this contamination by comparing the quark and 
antiquark contributions to the integral in Eq. fl41|) . We evaluated the antiquark contribu- 
tion using the same the cross section da q A N Y (x\) / dx\ that we found for an incident quark, 
Eq. (|45|), and using the simple expression 



0.7(1 -x) 8 

= ^ ( 69 ) 

for the atiquark distribution. For A = 184 and k = 2.28 GeV/fm, we find that the an- 
tiquark contribution is only 0.3% of the quark contribution at x\ = 0.1. Additionally, 
the relative size drops rapidly for larger values of x\, allowing us to conclude that the 
EMC/antishadowing effects make a tiny contribution to the DY cross section throughout 
the kinematical region we are considering and justifying our neglect of these effects in the 
data analysis. 

5.3 Why a previous analysis |T3|] found no energy loss 

One may wonder why the rate of energy loss —dE/dz deduced from the E866 data in a 
previous analysis Jl3|| is so much smaller (—dE/dz < OAAGeV/ fm) than the value we 



found fitting the same data. In this section we address this question. Q A value smaller than 
the string tension ~ 1 GeV/ fm would be at the same time both surprising and a serious 
problem for the theory. 

2 We compare only with the second of three scenarios for energy loss considered in |l3|. The first one 
(energy loss is proportional to the initial energy) violates the Landau-Pomeranchuk principle (see discussion 
in p"2|), while the third version (induced energy loss) corresponds to energy loss induced by multiple in- 
teractions of a quark propagating through infinite nuclear matter p6[. However the dominant contribution 



comes from vacuum energy loss (see Sect. 3.1) which has a constant and energy independent rate and is 



properly treated with the second version considered in fllJ 
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One obvious reason for a difference is the effect discussed in Sect. |3.2| and Fig. |6|, namely 
the observation that when the multiple scattering of the incident proton is properly consid- 
ered, the mean path of a quark available for vacuum energy loss in nuclear matter is more 
than twice as short as it would be if such scattering were neglected. Taking this effect alone, 
we estimate that it would increase the la upper limit for the rate of energy loss found in 
Ref. U] to -dE/dz < 1 GeV/fm. 



However, there are deeper reasons for the differences between the results of our analysis 



and that in Ref. |I3| , which treats the DY process as qq annihilation as in the leading-order 
parton model. These reasons are related to the fact that the parton model interpretation 
of the space-time development of the interaction is not Lorentz invariant and depends on 
the reference frame. As mentioned already in Sect. 2.1, the lack of Lorentz invariance arises 
because objects such as partons are quantum fluctuations rather than classical particles. 
Because of their quantum nature, the partons live and die differently in the different frames. 
Thus, one must avoid trying to look at the same parton in different reference frames, even 
though Xi is invariant under relative Lorentz boosts. For example, if the DY process is 
treated in the lab frame, then the fraction of light-cone momentum carried by the valence 
quark, x q , is larger than X\ [see Eq. (|9])]. However, in the dilepton rest frame x q — Xi, this 
is no surprise since they are not the same parton! 

5.3.1 Advantages of Target Rest Frame Formulation 

From the above discussion, it is clear that a consistent approach to the calculation of energy 
loss would, in principle, have to treat all effects in the same frame. In the case of the 
standard parton model approach jl], this would be the rest frame of the dilepton. In this 
frame, the moving nucleus is subject to a contraction by the Lorentz 7-factor, R* A = Ra/i-, 
where 

M , . 

7 = ~ • 70 

2 x 2 TUN 

The energy loss AE* in this reference frame is diminished by the same factor compared to 
the one in the nuclear rest frame, AE* = AE/j. Indeed, the vacuum energy loss is medium 
independent and proportional to L* = L/j. The induced energy loss is proportional to the 
Lorentz boosted nuclear density p* A = 7 pa, but also to (L*) 2 = L 2 /^ 2 . At the same time, 
the cloud of antiquarks in the target nucleon spreads in longitudinal direction over distance 
~ 2/M. Thus at small x 2 , 

(this is exactly the condition for long coherence length, i. e. for shadowing) the size of the 
antiquark cloud substantially exceeds the nuclear radius. 

Unfortunately, in the rest frame of the dilepton the position of qq annihilation cannot be 
localized to better than ~ 1/M, while the nucleus is squeezed by the Lorentz transformation 
to a pancake shape of a much smaller thickness ~ X2Ra^n I 'M . The annihilation process 
clearly lacks sufficient resolution to probe such small distances. Additionally, one cannot 
say whether it is the quark or the antiquark that propagates through the nucleus prior to 
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annihilation, a distinction that obviously must be made even before considering the issue of 
energy loss. 

For these reasons, it is quite tempting to switch reference frames (which, however, we 
have argued is invalid, unless it is done properly) and view the reaction from the rest frame 



of the nucleus. Just such a switch occurred in Refs. [11, 13|. For example, in Ref. [13 
the shift Ax\ caused by energy loss was first fit to the E866 data in the rest frame of the 
dilepton. The authors then switched to the rest frame of the nucleus, writing their Eq. (2) 



in this frame and assuming, in accordance with 0, [12), that k 2 is a constant independent 
of xi, x 2 and s. They find the rate of energy loss to be AE = AxxE/(3Ra/4), where 
E = 800GeV is the beam energy. 

One might try to avoid these inconsistencies by treating the DY reaction as qq annihi- 
lation in the rest frame of the nucleus, imagining the projectile parton traveling through 
nuclear matter and suddenly annihilating with a sea antiquark in one of the bound nu- 
cleons. This is, however, problematic for several reasons. First of all, the sea antiquarks 
available for such annihilation are defined only in the infinite momentum frame (this is why 
Weitzacker- Williams photons do not exist in the rest frame of the electron-they are its static 
field). Secondly, annihilation of a high-energy massless quark with an antiquark at rest to a 
dilepton of mass M violates energy-momentum conservation. To fix this problem, one may 
introduce next-to-leading order corrections, i.e. gluon radiation. As some of these gluons 
are radiated in advance of the dilepton and some later, one cannot treat the DY reaction 
as instantaneous but rather must consider its space-time development. Unavoidably, one 
arrives at the picture employed in present paper, where the DY pair is treated on the same 
footing as the gluons, i.e. as bremsstrahlung. 

One is thus led in a natural way to the target rest frame formulation, which clearly 
identifies the incident quark as the one that propagates through the nucleus and loses 
energy prior to the radiation of the 11 pair. The relative contribution of shadowing and this 
energy loss is governed by the coherence length appearing as an argument of the longitudinal 
form factor in Eq. (|64[). Specifying this dependence is essential for the determination of the 
correct value of the rate of energy loss from the data. This is because many of the E866 data 
are located at small x 2 , where the coherence length is relatively long (Fig.(|3|)), implying that 
shadowing dominates over energy loss in all but the heaviest nuclei. Ignoring this dependence 
on x 2 would thus tend to overemphasize shadowing and lead to a diminished rate of energy 



loss, just as in Ref. |13|| . 



5.3.2 Comparison of Results for Shadowing 

Still a substantial fraction of the E866 data is located at rather large x 2 where no shadowing 



is expected. For these events, another difference exists between Ref. [13] and our results. 



This can be traced to the fact that the shadowing correction in [13], which is taken from the 



phenomenological analysis by Eskola, et al. fbl, is quite different from ours in this region. 



We compare our shadowing calculated for tungsten at M = 4.5 depicted in Fig. [17| by the 
dashed curve, with the one from Ref. |T^] (solid thin curve) that was used in Ref. [TI| to cor- 
rect the E866 data. One can see quite clearly the difference between these two prescriptions. 
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Figure 17: Comparison of nuclear shadowing (dashed curve) and com- 
bined effect of shadowing and energy loss (thick solid curve) presented 
in our analysis, with "shadowing" for antiquarks found in jl^/ from the 
combined analysis of DIS and DY data (thin solid curve). All calcu- 
lations are done for tungsten and M = 4.5 GeV. The behavior of the 
dashed curve at small X2 is explained in the text. 



The source of the difference at large %2 arose, we believe, from a confusion between 
energy-loss effects and shadowing in the phenomenological extraction of the antiquark den- 
sities in Ref. |fL4 |. Shadowing for antiquarks cannot be extracted in a model independent 
way solely from DIS data, which are blind to the sign of the electric charge. This is why 



the E772 data for the DY reaction were included in the fit of Ref. |14| together with a 
variety of DIS data. Assuming QCD factorization, the evolution equations were applied to 
the parton distributions in nuclei. As usual, one needs to know the input x-distributions 
at a medium high scale Q$. The shapes of these distribution were guessed in Ref. Hl4 |, and 



then their magnitude was fit to data. No physical input beyond the QCD evolution was 
incorporated in the analysis. It is not a surprise that substantial "shadowing" effects were 
found for antiquarks even at large x 2 where no quantum interference is possible because 
of the shortness of the coherence length. It was assumed in Ref. [Il4l] that the ratio of the 



antiquark densities is constant for x > 0.08. The magnitude of this ratio was fit to the E772 
data and found to be less than 1, as one can see from Fig. 17. We think that this is where 
the misinterpretation of the energy loss effects as shadowing occurred in Ref. fTl 



as we emphasized in Sect. |5J] and one can see in Figs 
DY nuclear cross section even at small x\ (i.e. large x<i 



Indeed, 

IT] - [13], energy loss suppresses the 



where no shadowing is possible. 
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The thick solid curve in Fig. representing our result for the combined effect of energy loss 
and shadowing confirms this conjecture. Indeed, it is below the shadowing (dashed) curve 
at large £2 exactly by the same amount as the "shadowing" curve from Ref. |14|| . 

The influence of energy loss on the effective nuclear suppression of antiquarks found in 
14] extends down to smaller £2 (where most of the E866 data are concentrated), giving to the 
£2 distribution quite a different shape compared to the shadowing we calculated. In particu- 
lar, nuclear suppression is strongest at the kinematical limit of smallest £2 i n the E772/E866 
experiments, while shadowing shown by dashed curve vanishes in this limit. We should again 
emphasize that shadowing disappearing at small and large £2 is a model-independent ex- 
pectation since it is based only on the uncertainty principle and the kinematics of the DY 
reaction. 

Although the nuclear suppression from [I3| depicted in Fig. |I7] by the thin solid line 
exceeds the calculated shadowing effect shown by the dashed curve even at small £2 ~ 
0.015 — 0.02 this is not the reason for our disagreement with Ref. Indeed, as is stated 
above, we have checked that even with shadowing twice as strong as the one shown by 
the dashed curve in Fig. |T7| {i.e. stronger at small £2 than shadowing from Ref. |T^]), the 
resulting energy loss is essentially the same. We suppose that the more important effect is 
the different shape suggested in [|TJ]] for the £ 2 dependence of the nuclear suppression, which 
does not vanish at large £2. 

In conclusion, we believe that incorrect physical input at small £2 and the particular 

led to an incorrect result 



13 



antiquark density at large £2 employed in the analysis of Ref. 
of vanishing energy loss. Although this result was confirmed with much worse confidence 
by the recent analysis j|9| of data from the NA3 experiment, it comes as no surprise that 
they arrived at the same conclusions since their analysis was done essentially the same way 
as that in |13|j. 



6 Summary and Conclusions 

We have presented a new analysis of the nuclear dependence of the E772/E866 Drell-Yan 
lepton pair production data. This analysis makes use of a new formulation of DY in the 
rest frame of the target, according to which the lepton pair arises from the decay of a 
massive photon radiated by the incident quark in a bremsstrahlung process. We have been 
particularly interested in these data as a source of information on the rate of energy loss of 
a quark propagating through nuclear matter, encouraged by the observation that the rate 
of fall off of the DY cross section data with x\ near the point X\ = 1 is extremely rapid 
and therefore sensitive to the amount of energy lost by the quark prior to radiating the 
photon. Since shadowing has a similar effect to that of energy loss, namely suppression of 
the cross section at small £2, it is important that we have a good theoretical understanding 
of shadowing and that we explicitly use this understanding to model the shadowing. 

Identification of the dominant physical processes, as well as the determination of their 
relative contribution, is guided by the notion of coherence time, which can be interpreted as 
the lifetime of a virtual fluctuation of the incident quark into the massive photon (and final 
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quark). In particular, for short coherence times, the massive photon is released immediately 
from the fluctuation as the incident quark scatters from a bound nucleon. These processes 
clearly depend on the amount of energy lost by the soft incident quark as it propagates from 
its point of origin in a soft inelastic collision of the incident proton with a bound nucleon to 
the point at which it undergoes bremsstrahlung, and these processes are also sensitive to the 
energy lost by this quark for reasons given above. Because the coherence time is short, the 
photon is radiated immediately, and the final quark has no time to multiply scatter from 
other nucleons before the radiation process has come to completion. For long coherence 
times, the dominant processes correspond to terms in which the massive photon is radiated 
well in advance of any nuclear collision, and shadowing becomes an issue as the final quark 
begins to scatter from other nucleons of the nucleus. We describe shadowing quantitatively 
and without free parameters in terms of a phenomenological color-dipole cross section (fit 
to HERA data for the proton structure function), the nuclear thickness function, and the 
longitudinal nuclear form factor. Our model is shown to reproduce shadowing where it has 
been observed, namely in DIS on nuclei taken by the NMC collaboration. The transition 
between the regime of long and short coherence time is governed by the longitudinal form 
factor, with the longitudinal momentum transfer bearing a simple inverse relationship to 
the coherence length. 

Our theory, as described above, contains one unknown parameter, namely the rate of 
energy loss k of a quark propagating through nuclei. Using k as a free parameter, the theory 
was fit to the E772/E866 data for DY lepton pair production arising from bombardment of 
2 H, Be, C, Ca, Fe, and W targets by 800 GeV/c protons. The data were binned both in 
x\ and M, a very important issue for reliable discrimination between the effects of energy 
loss and shadowing. We found k = 2.28 ± 0.31GeV/ fm, in close accord with our theoreti- 
cal estimate of 2GeV/fm determined by considering string dynamics and gluon radiation. 
Numerous tests confirming the stability of our numerical analysis were made. They also 
provided a scale for the theoretical systematic uncertainty in the rate of energy loss, which 
we estimated at about ±0.5 GeV/fm. Since the value we found for k is substantially larger 
than the value reported in previous work, we pointed out the important differences between 
the two analyses. 

One should be cautious determining nuclear shadowing of sea quarks from DY data 
which may be substantially contaminated by energy loss. An essential demonstration our 
conclusions requires similar data for the DY reaction on nuclei at lower energies, where 
shadowing is of no importance but energy loss produces a stronger effect on the cross 
section. On the other hand, at much higher energies of RHIC and LHC one can completely 
disregard energy loss and test models for shadowing by direct comparison to data. We expect 
spectacular shadowing effects producing a strong suppression of the DY cross section for a 
wide range of x 1 . 
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